Formation of Rydberg macrodimers and their properties 
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We investigate the interaction between two rubidium atoms in highly excited Rydberg states, and 
find that very long-range potential wells exist. These wells are shown to support many bound states. 
We calculate the properties of the wells and bound levels, and show that their lifetimes are limited 
by that of the constituent Rydberg atoms. We also show how these /im-size bound states can be 
populated via photoassociation (PA), and how the signature of the ad- mixing of various ^-character 
producing the potential wells could be probed. We discuss how sharp variations of the PA rate 
could act as a switching mechanism with potential apphcation to quantum information processing. 

PACS numbers: 32.80.Rm, 03.67.Lx, 32.80.Pj, 34.20. Cf 
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In recent years, the strong interactions between Ryd- 
berg atoms, due to their exaggerated properties [H, have 
been detected experimentally and led to proposals 

for quantum computing e.g. to achieve fast quan- 
tum gates 0, Q or to study quantum random walks 0] ■ 
One effect of these interactions is the excitation block- 
ade 3, where a Rydberg atom prevents the excitation 
of nearby atoms [9|-[l3|. Recently, dipole blockade be- 
tween two atoms has been observed in microtraps 14 , 1^ , 
and a C-NOT gate implemented [31 . Another signature 
of these strong interactions is the molecular resonance 
features in excitation spectra, first observed in Rb 17 1 
and in Cs atoms [l8| . Other molecular features involving 
Rydberg excitations, such as the so called trilobite and 
butterfly states, where one atom_js in its ground state 
and another in a Rydberg state [l9|, have recently been 
detected [2^, while polyatomic molecules involving Ry- 
dberg electrons have also been studied [21, 
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In this article, we investigate the existence of long- 
range potential wells produced by two atoms excited into 
Rydberg states. As opposed to previous predictions of 
macrodimers [23] - doubly-excited Rydberg molecules - 
with very shallow wells due to induced van der Waals 
(vdW) interactions, the macrodimers we study here are 
due to the strong mixing of ^-characters of various Ry- 
dberg states. As we will see, these molecular wells are 
deeper (a few GHz), but still very extended with an equi- 
librium separation of 1 /im or so. We concentrate on 
wells existing in the vicinity of two rubidium (Rb) atoms 
excited to the np Rydberg state, for which the relevant 
molecular symmetries are 0^ , 0^, and 1„ [24| . 

We build our long-range molecular potential curves 
following the procedure described in |24| : we consider 
two free Rydberg atoms in states \a) = |n,£, j,mj) and 
\a') = \n' ,1' ,j' ,m'j), where n is the principal quantum 
number, £ the orbital angular momentum, and ruj is the 
projection of the total angular momentum j ~ £ + s 
onto a quantization axis (chosen in the z-direction for 
convenience) . The long-range molecular potential curves 
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FIG. 1. (Color online) 0^ molecular curves. Left panel: 
curves correlated to uili + n2l2 for n ~ 70, 69, and 68 
(top plot), 58 (middle) and 40 (bottom); black curves cor- 
respond to diagonalization near np -\- np asymptotes, and red 
curves near ns + ns. Right panel: enlargements near the 
70p3/2 + 70p3/2 asymptote (top: with a well correlated to 
69s -I- 72s), and the R dependence of lower asymptotes (mid- 
dle: n ~ 58; bottom: n ~ 40). Those last plots illustrate how 
the potential curves become flat as n decreases. 



in Hund's case (c) are calculated directly by diagonal- 
izing the interaction Hamiltonian consisting of both the 
Rydberg-Rydberg interaction and the atomic fine struc- 
ture 



24| . The basis states are built as follows: 

|a,a';r!g/„) ^ |a)i|a')2 - p{-lf+'' \a')^\a)2 



(1) 



where VL = nij + m'j is the projection of the total angu- 
lar momentum on the molecular axis and is conserved. 
The quantum number p, describing the symmetry prop- 
erty under inversion, is 1(-1) for g{u) states. For $7 = 0, 
the states' symmetry under reflection must also be con- 
sidered [13] . The basis ([T]) assumes no overlap of charge 
distributions between the two atoms, so the long-range 



2 



interaction can be expanded as an inverse power series of 
the nuclear separation distance R 24 1 . 

Fig. [1] shows the results of the diagonalized interaction 
matrix for the 0^ symmetry near various n£ + n£ asymp- 
totes. The left panel shows curves for several states lo- 
cated around np^^2 + '^^3/2 (in black) and ns + ns (in 
red) near n — 70, 69, 68, and lower asymptotes as well 
{n ~ 58 and 40). These plots depict the strong mixing of 
different ^-characters due mainly to the dipole-dipole in- 
teractions (although higher orders are also present). The 
right panel enlarges a major feature arising from the £- 
mixing near + 70^3/ 2, namely a large potential 

well (between and -2.5 GHz) with a depth ^ ^ GHz 
measured from its minimum (equilibrium separation) 
at roughly 40,500 oq (ao- Bohr radius), i.e. over 2 /im. 
The i?-dependence for lower asymptotes is also shown on 
the same panel; as n decreases the potential curves be- 
come flatter (for same R interval and energy scale). This 
is particularly true as n reaches 40 or so. 

The long-range well correlated to the (n — l)s + (ri + 2)s 
asymptote is a general feature for this system; we de- 
pict a few cases in Fig. HJa). The well depth De and 
equilibrium separation i?e scale roughly as n~'^ and v?'^ 
(Fig. mb)), respectively, in good agreement with the 
De oc and Re oc n^/^ scalings expected for a dom- 
inantly dipolar coupling between states. This last re- 
sult can be obtained by keeping only the major n- 
dependence, so that the interaction matrix between two 
asymptotes has the form hQ+n^ R~^ hi where /iq and 
hi contain the energy spacings between the asymptotes 
and the strengths of the dipole-dipole coupling, respec- 
tively. If we define p = R^^ and rewrite the inter- 
action matrix as n"^ [Hq + phi) , it becomes clear that 
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its eigenvalues Ei take the form Ei 



fi{p). Since 



dEi/dR = —3n'^R~'^fl{p) = at an extremum, we must 
have fl{p) =0, which is only fulfilled for some particu- 
lar values pe. In the parameter space of R and n, this 
leads to pe — n^R~^\R^ = const, so that Re oc n^/^, and 
fi{Pe) = n^Eili)^ = const, so that De oc n^"^. The differ- 
ence between the analytical and numerical scalings re- 
flects the more complex nature of the coupling than the 
assumed dipolar interaction. 

We return to the curve correlated to the 69s -I- 72s 
asymptote in the 0^ symmetry. Its well is deep enough 
to support several bound molecular levels, a few of which 
are listed in Table Uand shown in Fig. [2][c). Since the 
levels are separated by a few MHz, they could be detected 
by spectroscopic means. Oscillation periods of a few /xs, 
related to the MHz frequency of these deepest levels, are 
rapid enough to allow several oscillations during the life- 
time of the Rydberg atoms (roughly a few 100 /is for 
n ^ 70). The classical inner and outer turning points 
of the levels indicate very extended molecules, hence we 
keep the term macrodimer to describe them. 

The molecular electronic states \x\iR)) corresponding 
to the potential curves U\{R) are composed of several 
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FIG. 2. (Color online) 0^ wells correlated to the (n — l)s -|- 
{n + 2)s asymptote for n — 58, 70 and 85: (a) un-scaled, and 
(b) scaled to n = 58 with De ~ and Re ~ n^'^ . In (c), 
we show the deepest bound levels for n = 70 (original scale). 



molecular electronic basis states ([T]), the exact mixing 
(after diagonalization) varying with R: 



|XA(i?))=E#^(^)l^')- 



(2) 



Here, c^'^^(i?) are the eigenvectors and \j) the electronic 
basis states. As an example, we illustrate the compo- 
sition of the well correlated to the 69s -I- 72s asymp- 
tote in Fig. [31Ja). If we write \aa') = |a,a';0+), this 
well is composed primarily of |j) = |69s72s), |70s71s), 
|70p3/270pi/2), |69p3/271pi/2), and |69pi/271p3/2); these 
states are indicated by different colors in Fig.[3Ua). Using 
the same color labeling, we show the corresponding prob- 
abilities |cj(i?)p of all five states against R in Fig. [Sljb). 
As expected, the right side of the well (near the 69s -I- 72s 
asymptote) is composed primarily of the |69s72s) state, 
while the left side contains mostly |69pi/271p3/2), a state 
whose asymptote is above 69s -h 72s. Although one could 
have expected |68p72p) to provide most of the mixing, 
as suggested by Fig. [3]^ a) , the proximity of the 69s -I- 72s 
and 69pi/2 + 71p^/2 asymptotes and their strong dipole 
coupling result in this larger mixing. The probability for 
|68p3/272pi/2) is also included in Fig. |31^b). 

The coupling to states correlated to lower asymptotes 
could lead to predissociation of the macrodimers and 
their decay into free atoms, heating up the sample with 
more energetic collisions and ionization. To estimate 
the predissociation lifetime of our bound states, and be- 
cause of the large number of coupled electronic states, 
we assume that the total resonance width Fi associ- 
ated to the nonadiabatic transition from the molecular 
state 1 is just the sum over the widths Fi^ associated 
to transitions to individual adiabatic curves i, namely 
Fi = TiiTii. To compute the widths, we use a simple 
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TABLE I. Energies of the deepest bound levels (from the bot- 
tom of the well) and classical turning points for the well cor- 
related to the 69s + 72s asymptote of the 0^ symmetry. 
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Energy (MHz) 


Ri (a.u. 
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0.831242 


40,228 


40,679 
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2.499032 


40,068 


40,849 
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4.166823 


39,959 


40,970 
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5.824596 


39,870 


41,068 


4 


7.477361 


39,795 


41,154 


5 


9.125118 


39,728 


41,233 



two-channel approach in which the nonadiabatic cou- 
pling Vii between the close channel corresponding to 
the molecular bound state (with electronic adiabatic ba- 
sis state |xi(i?))), and the open channel corresponding 
to the dissociation state (with electronic adiabatic ba- 
sis state |x.(i?))), is Vu{R) = -^(Xi(i?)ljj Jj, 
where n is the reduced mass of the system. We then 
obtain the width Tu using a Green's function method 
Tii = 2Tr\{(/)y\Vu\(f)rcg)\'^ , where (/)rcg(i?) is the reg- 
ular, energy normalized solution of the open channel i 
(in the absence of channel coupling), and the bound vi- 
brational state 4>v{R) in the closed channel 1. For the 
69s72s Og potential curve, this approach give predissoci- 
ation lifetimes that are extremely long, at least 10^ years, 
and we thus can conclude that these long-range Rydberg 
molecules are limited by the lifetimes of the Rydberg 
atoms themselves. This result could be expected, since 
the most important states contributing to the well are 
above the 69s + 72s asymptote (Fig.[3l[b)). We note how- 
ever that for other wells due to mixing with lower asymp- 
totic electronic electronic states, predissociation would 
play a major role in the lifetime of macrodimers. 

In order to detect and probe the properties of the 
macrodimers predicted above, one can photoassociate 
two ground state atoms using intermediate Rydberg 
states that provide good overlap with the target long- 
range molecular states. Exciting the two ground state 
atoms into an intermediate Rydberg state ns, np, or nd 
will allow us to probe the different .^-characters (Fig. [3]) 
in the bound level. The photoassociation (PA) rate K„ 
into a bound level v can be calculated [2^ using 
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FIG. 3. (Color Online) (a) Potential well correlated to the 
69s + 72s asymptote, and (b) its composition, with probabil- 
ities |cj(i?)p of electronic basis \j) vs. R (inset: zoom of the 
inner region). Same color code used in (a) and (b): see text. 



R (corresponding to a flat curve), we can rewrite ([3]) as 

/>oo 2 

K. cx hh V |(dirf2)jf / dRcj):{R)c*{R)(f>g{R) , 

(4) 

where di = {nj£j\eri\ng£g) and d2 = {nj£j\er2\n'g£'g) are 
the electronic dipole moments between electronic states 
\ag,a'g) and \aj,a'j) for atom 1 and atom 2, respectively. 
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where h and I2 are the intensities of laser 1 and 2, 
\(j)v{R)) and \x\=i{R)) are the radial and electronic 
wave functions inside the well, respectively, \(f)g{R)) and 
\Xg{R)) are the radial and electronic wave functions of 
the ground state, respectively, and ri and e are the loca- 
tion and charge of the electron i. Using expression ^ 
for \x\=iiR)) , and assuming that \xg) is independent of 



FIG. 4. (Color online) PA rate vs. the detuning A from 
the 70p3/2 atomic state (a) on a linear scale, and (b) on a 
logarithmic scale. The intermediate state |41s, 41s) popu- 
lates the \np,n'p) character in the well (turquoise), and the 
|40p3/2, 4OP3/2) the |ns,n's) (red): (c) and (d) show the rate 
for the deepest levels (see text). 

We calculate K„ by assuming that the Rb 5s atoms 
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are first excited to an intermediate state that we con- 
sider to be the electronic "ground" states |xg), before 
being excited to their final electronic states. For sim- 
plicity, we consider lower asymptotes in the vicinity of 
n = 40 for which \xg) does not vary with R (Fig. [T]). As 
shown in Fig. [3J the final electronic state contain mostly 
\np^n'p) or \ns,n's); thus, we assume transitions from 
40p3/2 -I- 40p3/2 to the ns -\- n's components and from 
41s -f 41s to the np + n'p components. In Fig.|4j we show 
Kt, against the detuning A from the atomic 70p3/2 lev- 
els. For simplicity, we set (j)g{R) constant, although other 
choices could be made to enhance K^,, such as a gaussian 
corresponding to the motion of atoms trapped in optical 
tweezers. The PA rate starting from both atoms in 40p3/2 
is shown in red, and from 41s in turquoise, respectively, 
on a linear scale in (a) and a logarithmic scale in (b). 
In both cases, the rapid oscillation between a large rate 
for an even bound level {v = 0, 2, 4, . . .) and a small rate 
for odd levels {v = 1, 3, 5, . . .) gives the apparent envelop 
of the signal. This is illustrated in Figs. \^c) and (d), 
where a zoom of the deepest levels v is shown for both 
"ground" states (on a log-scale) . This can be understood 
from the integral in Eq.(j4]), which is near zero for an odd 
wave function <j)v{R)- 

The signature of a macrodimer would manifest itself 
by the appearance of a signal starting at A ~ —0.93 
GHz red-detuned from the 70p3/2 atomic level, and end- 
ing abruptly at ^ —1.36 GHz. In addition, the PA rate 
from either "ground" state 40p3/2 or 41s can reveal the 
details of the ^-mixing in the potential well. First, the 
overall larger signal for 41s reflects the scaling of c?id2 
appearing in Eq.(|3]), with 41s being closer than 40p3/2 to 
most atomic states appearing in the electronic molecular 
curve (see Fig. [3]). As A increases from -1.36 GHz, K„ 
mimic the probabilities |cj(i?)p. For 41s, the progres- 
sive decrease followed by sharp increases between —1.22 
and —1.16 GHz correspond to the slow decreases of most 
|np, n'p) components between i? ^ 40 — 50 a.u. and their 
sharp increases around i? ^ 33 — 35 a.u. (especially from 
|69pi/2, 7IP3/2), |69p3/2, 7IP1/2), and |70p3/2, 70pi/2); see 
Fig. [3]). This is followed by another large feature be- 
tween — 1.13 and —1.09 GHz (corresponding to the sharp 
increase in |69pi/2, 7IP3/2) and |69p3/2, 71pi/2) around 
R ^ 53,000 a.u.), and a "steady" rate before the sig- 
nal drops to basically zero. For 40p3/2, we obtain two 
obvious features: first, a significant drop in K„ between 
— 1.13 and —1.09 GHz which mirrors the decrease mainly 
in |69s, 72s), and second, a steady growth corresponding 
to the rising probability of |69s, 72s) beginning at about 
55, 000 a.u. (before the signal drops to basically zero). 

In conclusion, we have shown the existence of poten- 
tial wells for the 0+ symmetry of doubly-excited atoms 
due to ^-mixing. These wells support several bound lev- 
els separated by a few MHz with lifetimes limited by 
the Rydberg atoms themselves. These vibrational lev- 
els could be populated and detected by photoassociation. 



We showed that the signature of macrodimers would be 
the appearance of a signal beginning and ending abruptly, 
and that features in that signal could be used to probe the 
^-character of the potential well. The detection of such 
extended molecules with a lot of internal energy is in it- 
self a goal, and would allow to study how macrodimers 
are perturbed by surrounding ground state atoms. In 
addition to interesting exotic chemistry {e.g., by an ap- 
proaching Rydberg atom leading to a strong ultra-long 
van der Waals complex) , the dependence of the PA rate 
on the exact mixing (Fig.|l|) could potentially be used as 
a switch with application in quantum information {e.g., 
by turning on and off the blockade effect 

The work of N.S. was supported by the National Sci- 
ence Foundation, and the work of R.C. in part by the 
Department of Energy, Office of Basic Energy Sciences. 



[1] T.F. Gallagher, Rydberg Atoms (Cambridge University 

Press, Cambridge, 1994). 
[2] W.R. Anderson, J.R. Veale, T.F. Gallagher, Phys. Rev. 

Lett. 80, 249 (1998). 
[3] I. Mourachko et al, Phys. Rev. Lett. 80, 253 (1998). 
[4] M. Saffman, T. G. Walker, and K. M0lmer, 

larXiv:0909.4777v3 (2010). 
[5] D. Jaksch et al, Phys. Rev. Lett. 85, 2208 (2000). 
[6] I.E. Protsenko, G. Reymond, N. Schlosser, and P. Grang- 

ier, Phys. Rev. A 65, 052301 (2002). 
[7] R. Cote, A. Russell, E. E. Eyler, and P. L. Gould, N. J. 

Phys., 8, 156 (2006). 
[8] M.D. Lukin et al, Phys. Rev. Lett. 87, 037901 (2001). 
[9] D. Tong et al Phys. Rev. Lett. 93, 063001 (2004). 
[10] K. Singer ei a/.Phys. Rev. Lett. 93, 163001 (2004). 
[11] T. Vogt et a/.Phys. Rev. Lett. 97, 083003 (2006). 
[12] T. C. Liebisch, A. Reinhard, P. R. Berman, and G. 

Raithel, G. Phys. Rev. Lett 95, 253002 (2005). 
[13] R. Heidemann et a/.Phys. Rev. Lett. 100, 033601 (2008). 
[14] A. Gaetan et al Nature Phys. 5, 115 (2009). 
[15] E. Urban et al. Nature Phys. 5, 110 (2009). 
[16] L. Isenhower et al, Phys. Rev. Lett. 104, 010503 (2010). 
[17] S. M. Farooqi et al, Phys. Rev. Lett. 91, 183002 (2003). 
[18] K. R. Overstreet, A. Schwettmann, J. Tallant, and J. P. 

Shaffer, Phys. Rev. A 76, 011403 (2007). 
[19] C. H. Greene, A. S. Dickinson, and H. R. Sadeghpour, 

Phys. Rev. Lett. 85, 2458 (2000). 
[20] B. Vera et al. Nature 458, 1005 (2009). 
[21] I. C. H. Liu and J. M. Rost, Eur. Phys. J. D 40, 6571 
(2006); I. C. H. Liu, J. Stanojevic, and J. M. Rost, Phys. 
Rev. Lett 1 02, 173001 (2006); V. Bendkowsky et al, 
larXiv:09T2:4 058yl. 
[22] S. T. Rittenhouse and H. R. Sadeghpour, Phys. Rev. 

Lett. 104, 243002 (2010). 
[23] C. Boisseau, I. Simbotin, and R. Cote, Phys. Rev. Lett. 

88, 133004 (2002). 
[24] J. Stanojevic et al, Eur. Phys. J. D 40, 3-12 (2006); J. 

Stanojevic et al, Phys. Rev. A 78, 052709 (2008). 
[25] H. Friedrich, Theoretical Atomic Physics, Second Edi- 
tion, Springer- Verlag, 1998, ISBN 3540641246. 
[26] R. Cote and A. Dalgarno, Phys. Rev. A 58, 498 (1998). 



